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Based on the photometric properties of NGC 6822 we derive a robust star formation history. 
Adopting this history we compute 15 models of galactic chemical evolution. All of them match 
present-day photometric properties. The dark halo mass in all models evolves according to the 
mass assembly history predicted by a ACDM cosmology. We model the evolution of the baryonic 
mass aggregation history in this cosmological context assuming that part of the gas available for 
accretion never falls into the system due to two different physical processes. For seven models 
we assume that during accretion the universal baryon fraction is reduced by reionization only. 
The best model of this first group, a complex model with an early outflow, fits the observed 
gaseous mass, and the O/H, C/O, and Fe/O present-day values. This model requires a lower 
upper mass limit for the IMF than that of the solar vicinity, in agreement with recent results 
by other authors. We have also computed eight models where, in addition to reionization, the 
accreted baryon fraction is reduced by large-scale shock heating. The best model of this series, 
that also requires and early outflow and a lower upper mass limit for the IMF, can marginally fit 
the gaseous mass and the O/H and C/O observed values. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: individual NGC 6822 — galaxies: 
irregular — Local Group 
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ABSTRACT 



The main aim of this paper is to produce chem- 
ical evolution models of the irregular galaxy NGC 
6822. These models will be used to look for consis- 
tency with the cosmological context of galaxy for- 
mation as well as to study the possible presence 
of outflows from this object to the intergalactic 



In the field of chemical evolution often it is ar- 
gued that a model for a given galaxy is far from 
being unique. In this paper we use as many ob- 
servational constraints as possible to try to obtain 
a robust model with the idea of getting closer to 
a unique solution. The constraints used include: 
gas mass, chemical abundances, integrated colors, 
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and a star formation history derived from observa- 
tions. On the other hand the models used a gas in- 
fall history derived from the hierarchical standard 
ACDM cosmology. The models also used the IMF 
and the yields that reproduce the chemical evolu- 
tion of the solar neighborhood and the Galactic 
disk (Carigi et al. 2005). 

NGC 6822 is particularly suited for chemi- 
cal evolution models because its star formation 
history is well known (Wyder 2001, 2003), its 
present chemical composition might permit to de- 
cide if outflows to the intcrgalactic medium have 
been produced by this galaxy. NGC 6822, unlike 
the Magellanic Clouds, has not lost a significant 
amount of mass due to tidal effects. 

NGC 6822 is a member of the local group that 
apparently has not been affected by tidal effects 
from the Milky Way or M31 (Sawa & Fujimoto 
2005). There has been no recent interaction be- 
tween NGC 6822 and these two galaxies because: 
i) NGC 6822 is located 495 kpc away from our 
Galaxy and is receding from it at a radial veloc- 
ity of 44 km s- 1 (e.g. Trimble 2000), and ii) NGC 
6822 is also located at 880 kpc from M31, it is sep- 
arated from M31 by more than 90 degrees in the 
sky, and M31 is approaching our Galaxy at 121 
km s- 1 (Trimble 2000). 

An open problem in the study of the chemi- 
cal evolution of irregular galaxies is the cause of 
the difference between the observed effective oxy- 
gen yield and the oxygen yield predicted by closed 
box chemical evolution models. One solution to 
this problem might be due to the presence of O- 
rich galactic outflows (e.g. Garnett 2002). Sev- 
eral lines of reasoning indicate that O-rich out- 
flows produced by gas rich irregular galaxies are 
unlikely. Carigi et al. (1995, 1999), based on the 
C/O and O/H values for nearby irregular galaxies, 
have found that O-rich outflows have not played 
an important role in their evolution. A similar 
result has been obtained by Larsen et al. (2001), 
based on the N/O versus O/H relationship. These 
results indicate that O-rich outflows do not play 
a major role in the evolution of irregular galaxies. 
Another solution to the O yield problem might be 
due to the presence of dark matter (e.g. Carigi 
et al. 1999). 

The importance of outflows to the intergalac- 
tic medium from nearby irregular galaxies depends 
on many factors, like their total mass, the distri- 



bution in time and space of their star formation, 
and the presence of tidal effects (e.g. Legrand et 
al. 2001; Tcnorio-Tagle et al. 2003; Martin 2003; 
Fragile et al. 2004; Fujita et al. 2004, and refer- 
ences therein). 

In most galaxy chemical evolution models the 
infall of gas of primordial composition is modeled 
assuming a parametric form (e.g., Lanfranchi & 
Matteucci 2003). Here we will take a cosmolog- 
ical approach: the rate at which gas is accreted 
by the galaxy will depend on how the mass, dark 
and baryonic, falls into the galaxy (the mass as- 
sembly history, MAH), and the reionization his- 
tory of the universe. The MAH in turn depends 
on the assumed cosmology and the total mass of 
the galaxy (e.g., Avila- Reese et al. 1998; van den 
Bosch 2002). For a halo of a given mass, there is 
an ensamble of possible MAHs. In this paper we 
will take the mean MAH (defined in section 3.3) 
as representative of the MAH of NGC 6822. Be- 
fore the epoch of reionization, the history of the 
infalling gas is simply related by a constant to the 
MAH, afterwards the amount of gas accreted by 
the galaxy will depend on how fast the mass of 
the dark halo and the filtering mass grow with 
time (Gnedin 2000; Kravtsov, Gncdin, & Klypin 
2004). We will assume also that the gas infalling 
history is not altered by processes such as ram 
pressure or tidal striping. These assumptions are 
supported by the absence of warm-hot gas associ- 
ated to the Local Group and by the fact that this 
galaxy is isolated. 

In section 2 we define different masses present in 
an irregular galaxy and their role in chemical evo- 
lution models. The total mass and the halo mass 
assembly history within a ACDM cosmology are 
derived in sections 3.1 - 3.3. In section 3.4 we dis- 
cuss two scenarios for the baryonic mass building 
of the galaxy, one in which the reionization of the 
universe reduces the galaxy baryon fraction, below 
the universal one; and another in which in addi- 
tion to reionization the galaxy baryon fraction is 
also reduced by large-scale shock heating. In sec- 
tion 4 we estimate the gaseous mass and discuss 
the present day chemical composition to be fitted 
by the models. In section 5 we determine the vari- 
ation of the star formation rate and the luminosity 
with time. In section 6 we discuss the evidence in 
favor of gas outflows. In section 7 we present the 
general characteristics of the chemical evolution 
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models. In section 7.1 we discuss the large infall 
models, hereinafter L models, where the effect of 
gas shock heating is not considered. In section 7.2 
we discuss the small infall models, hereinafter S 
models, where the effect of gas shock heating is 
considered, this effect prevents a substantial frac- 
tion of baryons from reaching the galaxy. In sec- 
tions 8 and 9 we present the discussion and the 
conclusions. 

2. Definitions 

We can define the total mass, M to tai, as 

M tota i = M har + M DM , ( 1 ) 

where M bar is the baryonic mass, and M DM is the 
non-baryonic mass. 

Mbar can be expressed as 

M bar = M sub + M stars + M rem + M gas , (2) 

where M su b is the mass in substellar objects (with 
stellar mass, m, lower than 0.1 M©), and has also 
been called the baryonic dark matter, M star s is 
the mass of objects with m > 0.1 M , and M rem 
is the mass in compact stellar remnants. 
M gas can be approximated by 

M g as = M(H I) + M(He I) + M(H 2 ), (3) 

where we have neglected the ionized and the heavy 
element components. M gas is also given by 

Mgas = M 6ar (accreted) - M gas (outflow) - 

(M sub + M sta rs + M rem ) + Mr, (4) 

where Mbar (accreted) is the baryonic mass ac- 
creted from the intergalactic medium during the 
galaxy formation, M gas (outflow) is the baryonic 
mass ejected to the intergalactic medium during 
the galactic evolution, and Mr is the mass re- 
turned to the ISM by the evolution of the stars. 

We define the ratio of the baryonic mass ac- 
creted by the galaxy to the total mass as 

f acc = M &ar (accreted)/M tota; , (5) 



and the ratio of the baryonic mass that remains in 
the galaxy to the total mass as 

fgal = face ~ M gas (outflow)/M tota j . (6) 

In this paper, galaxy and halo are synonymous. 
3. Total Mass: M total 

3.1. Mtotai from direct observations of the 
rotation curve 

A high-resolution rotation curve that extends to 
5 kpc has been determined for this galaxy by Wel- 
drakc ct al. (2003) using the Australian Telescope 
Compact Array. Inside this radius, R, the mass of 
the galaxy given by Mi ower — is 3.5 x 10 9 M Q , 
where we are assuming that the circular velocity, 
V c , is equal to the rotational velocity determined 
observationally from the gaseous component of the 
disk and amounts to 55 km s _1 . With an observed 
baryonic mass of 4.3 x 10 s M0(see sections 4 and 
7) we find that NGC 6822 is a galaxy dominated 
by dark matter with Mbar/Mi ower — 0.12, inside 
R < 5 kpc. 

3.2. Mtotai from the cosmological context 

The total mass of a galaxy is a model dependent 
quantity. Here we will use a cosmologically moti- 
vated model to estimate the mass of the galaxy in 
question. Cosmological collisionless N-body simu- 
lations show that density profiles of dark matter 
(DM) halos can be described well by the NFW pro- 
file (Navarro, Frcnk, & White 1995, 1996, 1997) 

PNFw(r) = r/r s {llr/r s r (?) 

The characteristic radius r s is the radius where 
the logarithmic derivative of /9 NFW is —2, while 
the characteristic density po is given by p = 
4/0nfw(t's)- The NFW profile can also be char- 
acterized by a concentration parameter, c vir = 
Rvir/r s , and by the virial circular velocity V V i r — 
y/ G Mtotai /Rvir, where M to tai is the virial mass 
of the galaxy. The virial radius, R V i r , can be de- 
fined as the radius where the average halo density 
is S times the mean density of the universe ac- 
cording to the spherical top-hat model. Here 5 is 
a number that depends on epoch and cosmologi- 
cal parameters (fi ,^ A ); for a flat (0.3,0.7) ACDM 
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model 5 ~ 334 at the present epoch. In practice, 
we use cosmology to infer a relationship between 
the maximum circular velocity, V pea k, and V V i r , 
that is given by 



yl = 4 63 /(c OT r) 2 
v vtr r peak 



(8) 



(Bullock et al. 2001), where /(c OTr ) = ln(l + 
c v ir) — c v i r /(l+c v i r ). We thus can estimate the to- 
tal mass of the halo once V pea k and c V i r are given. 

We use as V pea k the rotational velocity at the 
outermost measured point. The mass of the halo 
may be underestimated because: i) the rotation 
curve for this galaxy may still rise beyond this 
radius and ii) we are neglecting any contribu- 
tion from the gas velocity dispersion. The con- 
centration of a DM halo depends on many pa- 
rameters such as nature of the dark matter, halo 
mass, epoch, normalization of the power spec- 
trum, etc. Here we will assume the "concordance" 
ACDM model (fi = 0.3, fl A = 0.7, h = 0.7, 
as = 0.9) where h is the Hubble constant in units 
of 100 km s _1 Mpc -1 and as is the rms of the 
mass fluctuations computed with the top-hat win- 
dow of radius 8 h^ 1 Mpc. If we take the median 
concentration estimated using the model by Bul- 
lock et al. (2001) and V pea k = 55 km s _1 from 
Weldrake et al. (2003) we find M total = 2.6 x 10 10 
M Q . This is the total mass that we will assign to 
NGC 6822 that implies an R vir — 76.3 kpc and 

Mi ower / Mtotal = 0.13. 

The adopted M to tai implies that baryons at 
present contribute to the galaxy mass budget with 
only 1.65% , such small values seem to be typical 
of dwarf galaxies. 

3.3. Mass assembly history: MAH 

Within the CDM galaxy formation paradigm a 
halo grows in a hierarchical fashion: self-bound 
structures form from bottom to top through the 
accretion and merging of smaller structures. A 
halo of a given mass today can arise from one 
of many possible mass aggregation histories. We 
use the recipe described in appendix A of van den 
Bosch (2002), which is based on the extended 
Press-Schechter formalism (Bond et al. 1991; 
Bower 1991; Lacey & Cole 1993), to derive the 
mean MAH for our given particular cosmological 
model and halo mass. 



van den Bosch (2002) finds, after experiment- 
ing with a variety of fitting functions, that the 
MAHs are well fitted by the following parametric 
form: 



log 



M(z) 



M 



-0.301 



log(l + z) 

log(l + Zf) 



(9) 



where M(z) and Mo are the mass of the halo at 
z and z = 0, respectively, and Zf and v are free 
parameters. According to this expression, Zf is 
defined as the redshift where half of the mass of the 
halo is acquired, and it is 1.50 for the total mass 
and cosmological context assumed in this paper. 
Based on Figure 6 of van den Bosch we took two 
additional models with extreme values of Zf, given 
by log(l + Zf) = 0.4 ± 0.1, to cover the scatter in 
the mass assembly histories. The mean MAH is 
computed assuming Zf — 1.50, the slow and fast 
MAHs are those computed adopting Zf = 1.00 and 
Zf = 2.16, respectively. Figure 1 shows the MAHs 
of these models (solid for mean and dotted lines 
for slow and fast). 
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Fig. 1. — The mean mass assembly history in units 
of the present-day mass (solid line) computed us- 
ing the recipe described in appendix A of van den 
Bosch (2002) with z f = 1.50 (the redshift where 
the halo acquires half of its mass) . The scatter in 
the diagram has been incorporated by using two 
extreme values for zj namely 1.00 and 2.16 (dotted 
lines). 
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3.4. M gas versus time 

A high fraction of baryons ~ 90% in the 
present-day universe do not reside in galaxies, 
the missing baryons (Fukugita & Peebles 2004). 
They can be in virialized regions of halos (hot gas 
in hydrostatic equilibrium) , in a diffuse phase with 
overdensities 8 < 1000 (5 = p/p — 1, where p is the 
mean density of the universe), or in a warm-hot 
phase with 10 5 < T(K) < 10 7 (Dave et al. 2001; 
Kang et al. 2005; Mo et al. 2005). This gas 
could have been expelled out of halos by strong 
winds or never became bound to halos in the first 
place. 

Photoionization of the IGM by quasars or early 
star formation is an efficient mechanism that im- 
pedes gas to fall into small galaxies. Below we 
estimate the effect that the reionization of the uni- 
verse could have had on the gas infalling history 
of NGC 6822. 

We use the formulas B1-B3 of appendix B 
of (Kravtsov, Gnedin, & Klypin 2004) to esti- 
mate the effect of reionization on the structure of 
NGC 6822. The formula B3 given below 



f B {M,z) 



ft 



[l + 0.26M F (z)/M] 3 ' 



(10) 



is an expression for the fraction of cold gas left 
in the halo of mass M (Gnedin 2000) and fb is 
the universal baryon fraction and amounts to 0.14. 
The filtering mass, Mp, defined as the mass of a 
halo which would retain 50% of its gas mass, is an 
increasing function of cosmic time and depends on 
when reionization starts, z r (see, for example, Fig. 
B10 of Kravtsov et al.). We use the mass assembly 
history M(z) of section 3.3 and equation (10) to 
compute the fraction of baryonic matter accreted 
by the galaxy at any given redshift, that is given 

by 



— I' 

M(z) J Zr 



M(z r ) 

M{z) 



fb 



+ 



f g (M,z)M'(z)dz, 



(11) 



where the prime denotes derivative with respect 
to redshift. 

In Figure 2 we plot f acc as a function of time for 
the mean MAH (solid line). The formula above is 
formally valid for z < z r , but it is clear that for 



redshifts greater than z r = 7 this fraction is just 
given by /{,. Once the epoch of reionization begins 
the amount of baryons in the galaxy at any time 
will depend on how Mp(z)/M evolves. Mp(z)/M 
varies by about a factor 2 from z = z r (t = 0.75 
Gyr) to z = (t = 13.5 Gyr, the present time) 
for the MAH model. Its rapid increase from z r to 
z ~ 4 (t <~ 5 Gyr) makes f acc decrease after z r . As 
time goes on, the drop in f acc is stopped because 
Mp(z)/M decreases. 

Yet, there is a second mechanism that might 
reduce the amount of baryons in galaxies: assume 
that a significant fraction of baryons was never in- 
corporated into the galaxy. Primordial gas could 
have been shock-heated by the collapse of non- 
linear structures of mild overdensities (pancakes 
and filaments), raising the entropy floor (e.g., Mo 
et al. 2005) and thus increasing the cooling time 
above the Hubble time. Although simulations 
with higher numerical resolution arc needed to 
confirm this prediction (Mo et al. 2005), present- 
day simulations seem to support this idea. Nu- 
merical simulations by Kang et al. (2005) show 
that about 40% of the gas at present is in a warm- 
hot phase (25% with 10 5 < T(K) < 10 7 and 15% 
with T(K) < 10 5 ). We will model this phenomena 
by simply assuming that this increase in tempera- 
ture in the intergalactic medium occurs suddenly, 
when the mass of the gas in the galaxy reaches 60% 
of its present-day value (in the absence of large- 
scale shock heating). From this epoch, z = 1.10 
(t = 5.37 Gyr), to the present the amount of ac- 
creted baryons in the galaxy does not increase. 
The second drop in f acc (see Figure 2, dotted line) 
is due to this effect. 

It should be clear by now that the evolution of 
face as well as its present-day value depend on the 
MAH. For the MAHs considered for NGC 6822 
and in the absence of large-scale shock heating, 
facc(z = 0) takes the following values: 0.033 (slow, 
z f = 1.00), 0.045 (mean, z f = 1.50), and 0.066 
(fast, Zf — 2.16). For completeness we also show 
in Figure 3 the evolution of the total (dashed line) 
and baryonic (solid line) mass for the mean MAH 
model, time is cosmic time. As in Figure 2, the 
dotted line represents the model where the large- 
scale shock heating of the gas is incorporated. 

At this time it may be convenient to summa- 
rize the steps we have taken to produce the gas 
mass accretion history. The infall of gas is propor- 
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Fig. 2. — Evolution of the ratio of the baryonic 
matter accreted relative to the total mass for the 
mean MAH, taking into account the effect of reion- 
ization (solid line), and in addition to reioniza- 
tion the effect of large-scale gas shock heating af- 
ter t = 5.37 Gyr (dotted line). Before t = 5.37 
Gyr the solid line represents both types of models 
since the effect of large scale shock heating occurs 
at t — 5.37 Gyr and does not affect the MAH for 
earlier times. We have assumed the onset of reion- 
ization at t r = 0.75 Gyr (z r = 7). Note that / acc 
equals the universal baryon fraction /& for t < t r 
and drops quickly below /& soon after and reaches 
a minimum. At t = 5.37 Gyr (z r — 1.10), which 
is the time when the mass gas reaches 60% of the 
total accreted mass of the large infall models, the 
curve divides in two: the solid and dotted lines. 



tional to that of the dark matter until the onset 
of the reionization epoch which we have taken to 
be z = 7. Once the reionization starts the accre- 
tion rate of gas is no longer proportional to that 
of the dark matter. In fact, in some models, mod- 
els with shock heating, the fall of gas is stopped 
while the accretion of dark matter continues. In 
this sense, the baryonic buildup of the galaxy 
is not certainly hierarchically. We will see later 
that the model favored by the observed chemical 
abundances requires an early galactic wind which 
in turn makes even less hierarchical the baryonic 
buildup of NGC 6822. This non hierarchical sce- 
nario has been called the downsizing problem and 
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Fig. 3. — Dark (dashed line) and baryonic (solid 
and dotted lines) mass evolution. The infall of gas 
follows that of the dark matter until the onset of 
the reionization of the universe. From then on, 
the gas mass growth is reduced in relation to that 
of the DM. After t = 5.37 Gyr in the small infall 
models, in which the intergalactic gas temperature 
increases due to large-scale shock heating (dotted 
line), the galaxy stops accreting gas once the gas 
mass reaches 60% of the total accreted mass of 
the large infall models (see main text); i. c., in 
the small infall models we assume that 40% of the 
gas will never be incorporated into the galaxy. 



has been discussed elsewhere (e.g. Arnouts et al. 
2005, Cowie et al. 1996). 

In section 7 we will describe two families of 
models for the increase of M gas versus time. Those 
models that accrete a relatively large M gas and 
consider the effects of reionization only, will be 
called large infall models (L models), and their 
chemical evolution will be discussed in section 7.1. 
Those models that accrete a relatively small M gas 
and that in addition to reionization consider the 
effect of shock heating will be called small infall 
models (S models) and will be discussed in section 
7.2. 
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4. Present Mass and Chemical Composi- 
tion of the Gaseous Component 

4.1. Chemical composition 

As observational constraints of the evolution 
models we will use the present chemical compo- 
sition of the interstellar medium of NGC 6822. To 
this effect we will combine the abundances of the 
supcrgiant stars and the H II regions. The super- 
giant stars were formed a few million years ago, 
therefore we expect them to have the same abun- 
dances than the H II regions. 

In Table 1 we present the adopted NGC 6822 
abundances. With the exception of the Fe abun- 
dance all the abundances for NGC 6822 are those 
derived by Peimbert et al. (2005) for Hubble V, the 
brightest H II region of the galaxy. The Fe abun- 
dance was derived by Venn et al. (2001) from 
two A supergiants and is in excellent agreement 
with the value derived by Muschielok et al. (1999) 
from three B supergiants. The gaseous Fe abun- 
dance derived from the H II regions can not be 
used as a constraint because most of the Fe atoms 
are trapped in dust grains. The O abundance is in 
very good agreement with the O /H values derived 
by Venn et al. (2001) for the two A supergiants 
mentioned before that amount to 12 + log O/H = 
8.36. 

There has been some discussion about the pos- 
sible presence of abundance gradients in irregu- 
lar galaxies, for example Pagel et al. (1978) state 
that a spatial abundance gradient, an analogous 
to those found in spiral galaxies, is small or ab- 
sent in the LMC and conspicuously absent in the 
SMC. The accuracy of the data on H II regions 
and supergiant stars for NGC 6822 is not yet 
good enough to establish the presence of a gra- 
dient. Works by Venn & Miller (2002) and Lee 
et al. (2006) are not conclusive, but seem to in- 
dicate that the O/H gradient is small or absent. 
NGC 6822 V is located about 1 kpc away from 
the galactic center and according to Gottesman & 
Weliachew (1977) about 70 % of the gaseous mass 
(M(HI) = 1 x 10 8 M©) is located within 2.6 kpc. 
In what follows we will consider that the chemical 
composition of NGC 6822 V is representative of 
the gaseous component of NGC 6822. 

Also in Table 1 we present the Orion and solar 
abundances for comparison. For Orion: Fe comes 



from B stars (Cunha & Lambert 1994) and the 
other elements from the H II region (Esteban et 
al. 2004). The solar abundances come from As- 
plund et al. (2005), with the exception of the He 
abundance that corresponds to the initial value 
computed by Christensen-Dalsgaard (1998) . 

The Orion nebula is more representative of the 
present chemical composition of the solar vicinity 
than the Sun for two main reasons: i) the Sun ap- 
pears deficient by roughly 0.1 dex in O, Si, Ca, Sc, 
Ti, Y, Cc, Nd and Eu, compared with its imme- 
diate neighbors with similar iron abundances (Al- 
lende Prieto et al. 2004), the probable reason for 
this difference is that the Sun is somewhat older 
than the comparison stars, and ii) all the chemical 
evolution models of the Galaxy predict a steady 
increase of the O /H ratio in the solar vicinity with 
time, for example the chemical evolution model of 
the solar vicinity presented by Carigi (2003), Ca- 
rigi et al. (2005), and Akerman et al. (2004) indi- 
cates that the O/H value in the solar vicinity has 
increased by 0.13 dex since the Sun was formed. 

4.2. Gaseous mass: M gas 

From the H I measurement by Huchtmeier & 
Richter (1986) and adopting a distance of 495 kpc 
we obtain an M(H I) = 1.36 x 10 8 M Q . From 
the helium abundance presented in Table 1 we ob- 
tain that M(Hc I)/M(H I) = 0.33. The M(H 2 ) 
estimated for NGC 6822 by Israel (1997) is about 
10 % of M(H I). Consequently, we will multiply 
the M(H I) + Af(Hc I) gaseous mass by a factor 
of 1.1 to take into account the contribution due 
to M(H 2 ). From the previous considerations we 
obtain that M gas = (1.98 ± 0.2) x 10 8 M . 

5. Luminosity Evolution and Star Forma- 
tion Rate 

5.1. Star formation rate: SFR and its nor- 
malization of the My 

The SFR used in chemical evolution models is 
given in general by a theoretical parametrization 
that depends on many physical properties that 
vary from object to object and in many cases are 
not well constrained. Consequently, in some cases 
the adopted SFR might not be realistic. Fortu- 
nately for NGC 6822 we have a SFR based on 
observations. 
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Table 1 

NGC 6822, Orion, and Solar Total Abundances 



Element 


NGC 6822 a 


Orion b 


Sun c 


12 + log O/H 


8.42 ±0.06 


8.73 ±0.03 


8.66 ±0.05 


log C/O 


-0.31 ±0.13 


-0.21 ±0.04 


-0.27 ±0.10 


log N/O 


-1.37 ±0.17 


-1.00 ±0.10 


-0.88 ±0.12 


log Fe/O 


-1.41 ±0.10 


-1.23 ±0.20 


-1.21 ±0.06 


log He/H 


10.91 ±0.01 


10.988 ± 0.003 


10.98 ±0.02 


Y 


0.2433 


0.2760 


0.2486 


Z 


0.0066 


0.0137 


0.0122 



a Muschielok et al. (1999); Venn ct al. (2001); Peimbert et al. 
(2005). 

b Cunha & Lambert (1994); Esteban et al. (2004). 
c Asplund et al. (2005); Christensen-Dalsgaard (1998). 



Wyder (2001, 2003) based on data obtained 
with the HST has derived the star formation rate 
{SFRi) of eight fields in NGC 6822. Since these 8 
fields correspond to regions of the galaxy with dif- 
ferent luminosities and taken together amount to 
20% of the total luminosity of the object (see bel- 
low), we will assume that they are representative 
of the whole galaxy during its evolution. Based on 
the SFRiS by Wyder we have constructed a SFR 
of the whole galaxy. 

We have to normalize the results by Wyder to 
the cosmological age adopted by us. Our cos- 
mological model indicates that the galaxy started 
forming 13.5 Gyr ago, we have assumed a SFR = 
for the first 1.2 Gyr, therefore the oldest stellar 
population has a representative age of 12.3 Gyr. 

Wyder (2001, 2003) only shows the SFH of 
region VIII for three assumed galactic ages (9, 12 
and 15 Gyr), for the other seven regions he only 
shows the SFH for an age of the galaxy of 15 Gyr 
(note that for the SFR time increases towards the 
future, while for the SFH time increases towards 
the past). In this paper we have assumed that 
the differences in the SFH between the ages of 
12 Gyr and 15 Gyr are small (as it is the case 
for region VIII), and have taken the sum of the 
SFH of all the regions for an age of 15 Gyr as 
representative of the whole galaxy. In addition we 
have adopted the SFH in the 15 to 7 Gyr age 



interval as representative of the 12.3 to 7 Gyr age 
interval. 

The total star formation rate of the eight fields 
will be called SFR(t) S f and will be given by 

8 

SFR(t) 8f =Y, SFR ^- ( 12 ) 

i=l 

With the SFR(t) 8f and the GALAXEV program 
of Bruzual and Chariot (2003) we obtain an My = 
— 14.28 at the present time (t = 13.5 Gyr), which 
implies that the eight HST fields produce 20% of 
the luminosity of the galaxy {My = —16.10 ± 0.1, 
Hodge 1977). We added to the GALAXEV pro- 
gram an additional subroutine that includes the 
metallicity evolution, Z(t) (Bruzual, private com- 
munication). The Z(t) used by us is given by the 
average of the [Fe/H] enrichment histories adopted 
by Wyder (2003 and private communication) un- 
der the assumption that [Fe/H] is proportional to 
the Z value. This Z(t) function will be called 

Zwyder{t)- 

The SFRtotai was computed from SFHs based 
on the Wyder code and was normalized based on 
the Bruzual-Charlot code. Both codes are dif- 
ferent, but we think that the procedure used to 
derive the luminosity evolution is reasonable, be- 
cause there are many common ingredients between 
GALAXEV and Wyder's codes. For example: 
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both codes use the same IMF (Salpeter 1955) in 
a similar mass range (from 0.1 M Q to 100 — 120 
M ), and the same stellar evolutionary isochrones 
(those of the Padova group dependent on metal- 
licity). 

Assuming that the rest of the galaxy behaves 
as the HST fields observed by Wyder we obtain 
that 



SFR total {t) = N SFR x SFR(t)i 



(13) 



where Nsfr is the normalization constant. For 
Nsfr = 5 we obtain My = —16.07 at the present 
time in excellent agreement with the observed 
value (My = — 16.10±0.1). In Figure 4 we present 
the SFRtotai obtained, this SFR to tai will be as- 
sumed in the models. 
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Fig. 4.— Star Formation Rate of NGC 6822 as a 
function of time determined from the star forma- 
tion histories of the 8 HST fields. 

The My was obtained with the GALAXEV 
code under the assumption of Salpeter IMF and 

Zwyder{t)- 

To reproduce the chemical evolution of the so- 
lar vicinity it has been found that the IMF by 
Kroupa, Tout & Gilmore (1993, KTG) is consid- 
erable better than the Salpeter IMF (e.g. Carigi 



et al. 2005). Therefore in the next section the 
chemical evolution models have been computed 
by adopting the KTG IMF. To test the robust- 
ness of the derived the SFR we have computed 
the My value with GALAXEV for the KTG IMF 
and three different Z(t) functions: Z\y y der(t), Z(t) 
of model 7L, and Z(t) of model 8S, our two best 
chemical evolution models (see section 7). For 
Zwyder{t), Z (t) of model 7L, and Z(t) of model 
8S, we obtained My = -16.22, My = -16.18, 
and My = —16.14, respectively; these values arc 
within one a of the observed My value, and im- 
ply that SFR presented in Figure 4 is robust in 
changes to the IMF and Z(t) functions. 

5.2. Luminosity and color evolution 

In Figure 5 and 6 we present the Mjj, Mb, and 
My absolute magnitudes evolution and the J — H, 
J — K, and V — J colors evolution, respectively, 
as a function of time for NGC 6822 given by the 
GALAXEV code with the KTG IMF and Z(t) of 
model 7L. 
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Fig. 5. — U, B, and V absolute magnitudes for 
NGC 6822 as a function of time. M v , M B , and 
My were computed from the SFR shown in Fig. 
4 using the GALAXEV code with metallicity evo- 
lution of model 7L (our best chemical evolution 
model) and the KTG IMF. 

In Table 2 we compare the observed magnitudes 
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and colors with those predicted by different models 
for the present time (13.5 Gyr). The magnitudes 
M B and My are taken from Wyder (2001), while 
the colors Mj - Mh and Mj - Mk from Jarret 
et al. (2003). The U - B color is from Hodge 
(1977) and the reddening correction for U — B was 
taken from Schlegcl ct al. (1998). The agreement 
between the model predictions and observations is 
excellent for the Mjj, Mb, and My magnitudes 
and the Mj - Mh and Mj - Mk color indexes 
(see also Figure 6). 
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Fig. 6. — Same as Fig 5, but for the J — H, J — K, 
and V — J colors. The observational constraints 
Mj - Mh and Mj - Mk are also presented for 
comparison (Jarret et al. 2003). 

It is not possible to compare the visual magni- 
tudes with the infrared magnitudes because the 
infrared observations did not include the whole 
galaxy (Jarrett et al. 2003). By adopting the 
My - Mj color of the Z{t) 7L model, we can 
determine the fraction of the integrated light of 
the galaxy observed in the infrared. The photo- 
metric evolution model predicts My - Mj = 1.32 
and Mj — —17.50. On the other hand the ob- 
served Mj value by Jarrett et al. (2003) amounts 
to —16.2 a value 1.3 mag fainter than predicted by 
the model. This difference implies that the aper- 
ture used by Jarrett et al. (2003) included only 30 
% of the integrated light of the galaxy, in agree- 
ment with the comment in their section 6.5 where 



they mention that in general they observed only 
from 20 % to 50 % of the integrated light of dwarf 
irregular galaxies. 

From Table 2 it can be seen that the different 
models agree with the observations at one o level. 
In other words the absolute luminosity in different 
bands as well as the color indexes are almost in- 
dependent of the three different Z{t) functions as 
well as the two IMFs used and consequently that 
the SFR presented in Fig. 4 is robust and can be 
used for all the chemical evolution models. That is 
the SFR has been derived from observations and 
it will be kept fixed for all the chemical evolution 
models. 

6. Galactic Winds 

A common hypothesis in chemical evolution 
models of dwarf galaxies is to assume that these 
galaxies expelled part of their gas during their life- 
time: a galactic wind is mostly introduced in the 
modelling to explain the low chemical abundances 
found in these type of galaxies and NGC 6822 
seems to be no exception. From an observational 
point of view the amount of evidence of outflows in 
star forming galaxies is increasing with time and 
we refer the reader to Veilleux et al. (2005) for 
an updated discussion of galactic winds. On the 
other hand, from a theoretical point of view, nu- 
merical simulations show that a blow-away outflow 
situation is not rather uncommon phenomenon in 
low-mass galaxies (e.g Fujita et al. 2004). 

In order to see if NGC 6822 could have had a 
galactic wind during its lifetime, specially at high 
redshift, we have computed its present date bind- 
ing energy as well as the thermal energy history 
of the gas. In Figure 7 we present the thermal en- 
ergy of the gas as a function of time assuming no 
outflows. 

The thermal energy was computed as follows: 
we define Ether as the thermal energy of the gas 
from supernovae, 



ETHER{t) 



-f 

Jo 



e(t - T)RSN(T)dT (14) 



where RSN is the supernovae rate and includes 
types la and II SNe. The e(t) function represents 
the evolution of the thermal energy in the hot, 
dilute interior of the supernovae remnants and is 
given by 



10 



Table 2 

Present-day Photometric Properties for different Initial Mass Fuctions and Metallicity 

Evolutions 



IMF 


Salpeter 


KTG 


KTG 


KTG 


Observations 


Z(t) 


Wyder 


Wyder 


Model 7L 


Model 8S 




Mu 


-15.97 


-15.91 


-15.82 


-15.74 


-15.83 ± 0.16 


M B 


-15.74 


-15.81 


-15.76 


-15.70 


-15.60 ± 0.13 


M v 


-16.07 


-16.22 


-16.18 


-16.14 


-16.10 ± 0.10 


Mj - M H 


0.52 


0.53 


0.55 


0.56 


0.49 ± 0.04 


Mj - M K 


0.66 


0.66 


0.70 


0.72 


0.63 ± 0.05 


My - Mj 


1.18 


1.26 


1.32 


1.35 




Znow 


0.0024 


0.0024 


0.0069 


0.0084 


0.0066 ± 0.0011 



it) 



0.72e 

0.22e (t/tc)-° M2 



if t < t c , 
if t > t c , 



where eo = 10 51 erg and tc is a cooling time scale 
taken from Cox (1972). For further details see 
Carigi, Hernandez & Gilmore (2002). 

Also in Fig. 7 we show the binding energy of 
NGC 6822 for the present gas content based on 
the following equation: 



E 



BIND 



G 



M gas M tota i 

Ryi r 



(15) 



As can be seen from Figure 7 the thermal en- 
ergy available, released during the lifetime of the 
galaxy, is higher than the present binding energy 
of the gaseous content of the galaxy, which im- 
plies that part of the thermal energy was lost by 
the galaxy in the past. This result is in agreement 
with the chemical evolution models that require 
a substantial gas outflow to match the present 
gaseous content with the observed O /H value (see 
section 7). 

As we mention above, the computation of the 
thermal energy content of NGC 6822 is based on 
the cooling scheme by Cox (1972). A more elab- 
orated treatment of the SN remnant cooling has 
been presented by Tantalo ct al. (1998). They 
computed the cooling time as a function of gas 
metallicity and found that for metallicities below 
solar it is larger than that predicted by Cox. Thus, 



in absence of winds, the thermal energy in NGC 
6822 should lie above the one shown in Fig. 7. 

The outflows needed to explain the excess ther- 
mal energy predicted by the observed SFR and 
by the chemical evolution model presented in sec- 
tion 6 very likely occurred in the first few Gyr 
of evolution (t < 5 Gyr) because: i) To propiti- 
ate an efficient blow away process strong bursts 
of star formation are needed (Tenorio-Tagle et al. 
2003), from the SFR it is not possible to rule out 
the presence of several outburst processes at early 
times (t < 5 Gyr), because the SFR has been av- 
eraged at long time intervals (~ 2.5 Gyr), and it is 
possible for one or two large bursts of star forma- 
tion, to go unnoticed in the derived SFR, on the 
other hand, at later times the SFR has been av- 
eraged for shorter time intervals (~ 0.5 Gyr), and 
it is more difficult for a high SFR to go unno- 
ticed; ii) in the past, in a hierarchical cosmology, 
the potential wells of galaxies, for a given present- 
day mass, are shallower; iii) strong star formation 
events are expected to occur because of the higher 
galaxy merger rate (e.g., Mailer et al. 2005); iv) 
after an outflow occurs, M gas gets diminished and 
the possibility of a second outflow due to a burst 
of star formation increases because there is a lower 
amount of gas in the ISM to reduce the velocity 
of the SN ejecta, and v) the low metallicity at 
early times increases the cooling time relative to 
the classical Cox (1972) value, and therefore in- 
creases the thermal energy available propitiating 
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Fig. 7. — Gaseous thermal energy predicted by 
SFR under the assumption of no gaseous out- 
flow from the galaxy based on the formulation by 
Cox (1972). We also show the binding energy of 
the present gas content of the galaxy. The differ- 
ence between the predicted thermal energy and the 
binding energy at present implies that the galaxy 
must have lost thermal energy due to outflows in 
the past. 

the advent of galactic winds. As can be seen from 
our chemical evolution models, for example mod- 
els 7L and 8S (see section 7) the Z value for t < 5 
Gyr is about 0.001, thirteen times smaller than 
solar. 

It is beyond the scope of this paper to study 
the outflow rate as a function of time. 

7. Chemical Evolution Models 

We have computed chemical evolution models 
with a very low number of free parameters. The 
cosmological context provides the accretion rate, 
and the luminosity model provides the SFR as a 
function of time, both in an independent way. The 
remaining ingredients related with stellar popula- 
tions: the IMF, the stellar yields, and the per- 
centage of binary systems which explode as SNIa, 
are inferred from a successful chemical evolution 
model of the solar vicinity and the Galactic disk 
by Carigi et al. (2005). 



To construct the models for NGC 6822, we have 
relied upon the model of the solar vicinity by Ca- 
rigi et al. (2005), consequently we have adopted 
the following assumptions: i) the IMF proposed 
by Kroupa, Tout, & Gilmorc (1993, KTG) for 
the 0.01 - 80 M© range, ii) for massive stars, 
those with m > 8 M©, we have used yields that 
depend on the initial metallicity as follows, for 
Z t = yields by Chieffi & Limongi (2002), for 
Z t = 1 x 10~ 5 and 0.004 yields by Meynet & 
Maeder (2002), and for Z t = 0.02 yields by Maeder 
(1992), iii) for low-and-intermediate-mass stars, 
those with m < 8 M© and from Zi = 0.004 to 
Zi = 0.02 we have used the yields by Marigo, 
Bressan, & Chiosi (1996, 1998), and Portinari et 
al. (1998), iv) for Type la supernovae, we have 
used the yields for Zi = 0.02 by Thielemann et 
al. (1993), these yields are almost independent 
of metallicity, and v) the fraction of binary sys- 
tems with total mass between 3 < m(M©) < 16 
that becomes SNIa (Abin) amounts to 0.05. For 
each set of yields, linear interpolations for differ- 
ent stellar masses and metallicities were made. For 
metallicities higher or lower than those available 
we adopted the yields predicted by the highest or 
lowest Z available, respectively. 

In some models, we will adopt two different 
types of outflows: i) we call "selective outflows" 
those constituted only by material ejected by su- 
pernovae of Type II and Type la, in this case 
all the material ejected by massive stars (Type 
II SNe) is lost to the IGM, and ii) we call "well 
mixed outflows" those that reach the IGM after 
the ejecta by MS, LIMS, and Type la SNe have 
been completely mixed with the ambient ISM. We 
have assumed that the rate of the well mixed out- 
flow is proportional to the SFR and that the pro- 
portionality coefficient is constant during the pres- 
ence of the outflow. With both types of outflows 
the galaxy loses material that never gets back. 
The duration of the outflow is another free pa- 
rameter. 

We computed a large set of chemical evolu- 
tion models. From these, we will discuss in de- 
tail 15 models in sections 7.1 and 7.2. All models 
evolve with the mean MAH computed in section 
3.3. Seven of the 15 models follow a baryonic mass 
accretion history that does not assume large-scale 
heating of the gas, while the other eight models do. 
In some models we included outflows of baryonic 
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matter. All models fit the observed colors evolu- 
tion in Section 5.2, and twelve of them have the 
same values for M su b, M stars , and M rem , at the 
present time, these values amount to 3.48 x 10 7 
M , 1.73 x 10 s Mq, and 2.31 x 10 7 M Q . The 
exceptions are Model 7L, 7S, and 8S, that have 
different upper mass limits for the IMF, for these 
models the values of M su b, M sta rs, and M rem are 
practically the same values than those predicted 
by the other models and they will be presented in 
sections 6.1 and 6.2. The Mb ar determined from 
the observed gaseous mass and the other terms 
present in eq. 2 amounts to 4.3 x 10 8 M . 

7.1. Large infall models 

We have computed several models that adopt 
all the ingredients presented in the previous sec- 
tions, specially in Sections 3, 5.1 and 7. In par- 
ticular, in Models 1L - 7L we consider that the 
baryonic mass accreted by NGC 6822 is that sug- 
gested by the cosmological context without large- 
scale shock heating. 

In Table 3 we present the main characteristics 
of seven models. Columns 2 and 3 describe the 
type of outflow and its duration. Columns 4 to 
8 present five properties predicted by the models: 
M gasi O/H, C/O, N/O, and Fe/O, as well as their 
observed values. The M gas and the O/H values 
are the most critical observational constraints for 
NGC 6822, therefore a model that fits the C/O, 
N/O and Fe/O observed values, but that does not 
fit the observed O/H value, is a poor model be- 
cause it does not fit the C/H, N/H, and Fe/H ob- 
served values. 

In Model 1L we assume that all the baryons 
available after the reionization fall into the galaxy 
and remain in the galaxy during its whole evolu- 
tion. In Figure 8 we show the evolution history of 
six properties of the model up to now. Model 1L 
grossly fails to fit M gas and O/H, since it predicts 
4.7 times more gas than observed and a lower O/H 
value than observed by 5 a (see Table 3 and Figure 
8). 

In order to reduce the enormous difference be- 
tween the predicted M gas value by Model 1L and 
the observed value we have computed models with 
well mixed outflows during the early history of the 
galaxy. 

Model 2L is built to reproduce the M gas ob- 
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Fig. 8. — Model 1L. Mass formation history and 
predicted chemical evolution. This model assumes 
that all the gas available after the reionization 
falls into the galaxy and does not include out- 
flows. The present-day observational constraints 
(filled circles) are also presented for comparison. 
The model fails drastically. 

served value under the assumption of an instanta- 
neous well mixed outflow. This event had to occur 
at 6.3 Gyr, because at that time the M gas of the 
galaxy amounts to the difference between the ob- 
served and the predicted M gas by Model 1L (7.38 
xlO 8 M Q ). This model reproduces M gaSl C/O 
and Fe/O, but the O/H and N/O values are higher 
than observed by 2.0 a and 1.2 a respectively (see 
Table 3). O/H is higher because according to the 
SFR most of the stars are formed after the out- 
flow event and the O ejected by the stars to the 
ISM is kept by the galaxy 

Based on the models by Fujita et al. (2004) we 
find that in order to produce an instantaneous well 
mixed wind in NGC 6822 at 6.3 Gyr, when its to- 
tal mass is 1.7 x 10 10 M Q , about 10% of the gas 
should be converted into stars. This amount of 
gas corresponds to ~ 6 x 10 7 M Q which accord- 
ing to the derived SFR is 2.5 times higher than 
the mass converted into stars by NGC 6822 in an 
interval of 1 Gyr centered at 5.8 Gyr (see Figure 
4) . Consequently, we consider Model 2L to be un- 
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likely. 

Model 3L is built to reproduce the O/H ob- 
served value under the assumption of a well mixed 
outflow that starts at the same time than the star 
formation (t = 1.2 Gyr). To reproduce the ob- 
served O/H value the model needs the outflow 
phase to last 4.0 Gyr. During the outflow phase 
the mass lost to the IGM is the maximum possi- 
ble, that means that the gas left is the minimum 
required to maintain the SFR given in Fig. 4. 

Model 3L reproduces all the observational con- 
straints related with chemical abundances, but not 
the M gas observed value: the M gas predicted is 
higher than observed, by a factor of 1.65 (8 a) 
(see Table 3 and Figure 9). The gas mass lost is 
6.09 xlO 8 M Q , 8.5 times the mass of stars formed 
during the first 4.0 Gyr. During the wind phase 
the outflow rate is equal to 7.51 times the SFR. If 
it turns out that the observed gas mass value has 
been underestimated (which is unlikely since most 
of the gas is in the form of HI that is observed), 
this model might become relevant. 
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Fig. 9. — Model 3L. Mass formation history and 
predicted chemical evolution. This model assumes 
the gas mass assembly history of the large infall 
models and includes a well mixed outflow during 
the first 4 Gyr of the star formation. 

Model 4L is built to reproduce the M gas ob- 
served value under the assumption of a well mixed 



outflow. This model requires to lose more gas than 
Model 3L, therefore the duration of the outflow 
has to be longer. To fit M gas we need the out- 
flow to last 5.1 Gyr. During the outflow phase 
the mass lost to the IGM is the maximum possi- 
ble, that means that the gas left is the minimum 
required to maintain the SFR. During the wind 
phase the outflow rate is equal to 7.11 times the 
SFR. The gas mass lost is 7.35 xlO 8 M Q , 8.0 
times the mass of stars formed during the first 
5.1 Gyr of the star formation. The predictions of 
this model are identical to those of Model 2L (see 
Table 3), from ~ 7 Gyr up to now, because: i) 
both models at ~ 5 Gyr had lost practically all of 
the gas and started accumulating gas again at the 
same rate, and ii) both models had the the same 
SFR after 5 Gyr. In other words, after the end 
of the outflow phase both models have the same 
stellar populations that pollute equally the same 
amount of gas. 

Model 4L agrees with the observed M gas value 
but overproduces O, therefore we will modify 
Model 4L to fit the O/H ratio and the other 
observational constraints by adopting extra as- 
sumptions in Models 5L, 6L, and 7L. 

In Model 5L we have assumed that 0.5 Gyr ago 
a primordial gaseous cloud of 0.34 xlO 8 M was 
accreted by NGC 6822. We have assumed the least 
massive cloud needed to reproduce the O/H value. 
With this extra infall, the model can reproduce all 
the observational constraints related to chemical 
abundances with the exception of N/O, but the 
predicted M gas is higher by a factor of 1.2 (see 
Table 3). The accretion of a huge gaseous cloud of 
primordial material recently is more unlikely than 
at earlier times. If we assume the accretion of a 
primordial cloud at earlier times (t < 13 Gyr), 
the model would require a more massive cloud, in- 
creasing the disagreement with the observed M gas . 

To assume the infall of a cloud of 3.4 x 10 7 M 
of primordial gas at t = 13 Gyr is unlikely for the 
following reasons: i) this gaseous cloud has about 8 
% of the baryonic mass of the galaxy and therefore 
could have formed part of the NGC 6822 satellite 
system, ii) at such late time a gas cloud of this 
mass could have had time to form stars and these 
stars could have changed the primordial chemi- 
cal composition of the gas cloud (according to De 
Block & Walter (2005) the disk of NGC 6822 does 
have an embedded gaseous cloud with comparable 
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mass but with an old stellar population), and iii) 
the total gaseous mass accreted by the galaxy in 
the last 0.5 Gyr according to the MAH amounts 
only to 5 x 10 5 M Q . Consequently, we consider 
Model 5L unlikely. 

In Model 6L we have assumed a second outflow 
of the selective type, during the last 1 Gyr, when 
the SFR is most important. In this outflow, all 
the material ejected by SNII and SNIa is lost to 
the IGM. If the SNIa ejecta remain in the galaxy, 
the predicted Fe/O value becomes —1.30 dex, in 
marginally agreement with observations (see Table 
3). Model 6L is unlikely because between 12.5 and 
13.5 Gyr the SFR does not indicate the presence 
of a major starburst (see Figure 4) needed to drive 
the outflow, and the potential gravitational well is 
largest. Moreover, we think that it is unlikely that 
the material ejected by every SNe in the last Gyr 
escapes to the IGM. 

In Model 7L we have reduced the IMF mass up- 
per limit, m up , from 80 M Q to 60 M Q . When m up 
is reduced, the percentage of MS, main produc- 
ers of O, is decreased by a factor of 1.18, and the 
predicted O/H value is reduced by 0.12 dex. This 
fact produces an increase of the N/O and Fe/O 
predicted values to -1.06 dex and -1.24 dex, re- 
spectively, more than 1 a away from the observed 
values. The Fe/O value can be adjusted by as- 
suming a lower fraction of binary systems that be- 
come SNIa (Abin). For this reason we have re- 
duced Abin from 0.05 (the value assumed in the 
rest of the L models) to 0.02. Model 7L fits all the 
observational constraints with the exception of the 
N/O ratio (see Figure 10 and Table 3). 

With m up = 60 M Q , and Abin = 0.02, the stel- 
lar populations arc not the same than those in 
the solar vicinity. With these values the model 
by Akerman et al. (2005) or Carigi et al. (2005) 
cannot reproduce the chemical observational con- 
straints for the solar vicinity, implying that the 
IMFs of the solar vicinity and of NGC 6822 are 
not the same. Other authors have also argued in 
favor of smaller m up values for dwarf galaxies, for 
example, Weidner & Kroupa (2005) suggest that 
m up is lower in dwarf galaxies than in the solar 
vicinity; the smaller m up implies that the num- 
ber of Type II SN per generation of stars has to 
be smaller than in the solar vicinity. Their re- 
sults depend on the SFR. For a continuous SFR 
during the last 14 Gyr, their IMF slope is steeper 



5 10 15 7.5 8 8.5 




Fig. 10. — Model 7L. Mass formation history and 
predicted chemical evolution. This model assumes 
the gas mass assembly history of the large infall 
models and includes a well mixed outflow during 
the initial 5.1 Gyr of the star formation. The IMF 
m up = 60 M Q , and the fraction of binary systems 
that becomes SNIa is 2 %. 

than the Salpeter IMF for massive stars, as in the 
case of the IMF assumed in this paper, and has 
an m up ~ 10 M Q . For a single 100 Myr burst 
their slope is like the Salpeter's one and has an 
m up <~ 100 M . The SFR computed by us is 
nearly constant during the first 5 Gyr and then 
shows several bursts. Based on the work by Wei- 
dner & Kroupa the IMF adopted by us should be 
dependent on time, but that is out of the scope 
of this paper. We consider that a constant IMF 
with m up — 60 M© is a good approximation to 
the IMF proposed by Weidner & Kroupa (2005). 
Other authors also argue that galaxies with a lower 
SFR have a lower m up value in agreement with 
our suggestion (e.g. Goodwin & Pagel 2005, and 
references therein). 

The change in m up from 80 M© to 60 M Q af- 
fects the O/H evolution, but it docs not affect the 
photometric evolution nor the thermal evolution. 
We have computed a photometric evolution model 
with m up = 60 M© and obtained Mjj = —15.82, 
M B = -15.77, M v = -16.19, Mj - M H = 0.55, 
Mj - M K = 0.70, and M v - Mj = 1.32 practi- 
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Table 3 

Inputs and outputs of large infall models 



Model 


Wind 




M gas 




Abundances 






Type Duration 


(108 M ) 


12+log(0/H) 


log(C/0) 


log(N/0) 


log(Fe/0) 






(Gyr) 












1L 






9.36 


8.13 


-0.35 


-1.22 


-1.43 


2L 


Well Mixed 


0.1 


1.99 


8.54 


-0.33 


-1.16 


-1.35 


3L 


Well Mixed 


4.0 


3.27 


8.42 


-0.35 


-1.19 


-1.38 


4L 


Well Mixed 


5.1 


2.01 


8.54 


-0.33 


-1.16 


-1.35 


5L a 


Well Mixed 


5.1 


2.34 


8.48 


-0.34 


-1.17 


-1.35 


6L b 


Well Mixed 


5.1 


1.98 


8.42 


-0.27 


-1.09 


-1.37 


7L c 


Well Mixed 


5.1 


2.01 


8.42 


-0.24 


-1.06 


-1.41 


Obs. 






1.98 ± 0.20 


8.42 ± 0.06 


-0.31 ± 0.12 


-1.37 ± 0.18 


-1.41 ± 0.10 



a Extra 3.4 x 107 M of primordial gas accreted at 13.0 Gyr. 

b All the material expelled by Type II and la SNe during the last Gyr is ejected to the IGM. 

c The upper mass limit for the IMF amounts to 60 M Q . The fraction of binary systems that become SNIa 
amounts to 0.02. 



cally identical to those derived with m up = 80 M 
presented in Table 2, therefore the change in the 
SFR due to the reduction of from 80 M Q to 
60 Mq is negligible. Similarly, the gaseous thermal 
energy produced by the SFR under the assump- 
tion of no outflows diminishes from 1.65 xlO 55 erg 
to 1.63 xlO 55 erg when m up changes from 80 to 
60 M . 

For Model 7L the current values of M suo , 
M s tars, and M rem , amount to 3.49 x 10 7 M Q , 
1.73 x 10 8 M , and 2.31 x 10 7 M©, respectively. 

It can be shown that the use of the Salpeter 
IMF would imply a smaller m up than that required 
by KTG IMF to fit the observed O/H value. Con- 
sequently since all of the O is produced by massive 
stars and only a fraction of the C is produced by 
them, the C/O value would become considerably 
larger than observed (e.g. see model 7S). 

In Table 4 we present two carbon budgets: 
the stellar production one, and the interstellar 
medium one for L models. The production is 
due to massive stars (MS), low and intermediate 
mass stars (LIMS), and SNIa. Also in Table 4 
we present for comparison the carbon budgets for 
the solar vicinity (Carigi et al. 2005). For Model 
1L both budgets agree because there is no outflow 
present, a similar situation prevails for the solar 
vicinity. For Models 4L, 6L, and 7L the contribu- 
tion of massive stars to the C abundances in the 
ISM is smaller than the stellar production due to 
outflows to the IGM. Since for Models 4L and 6L 
the IMF and the SFR are the same, the higher 
contribution of MS to the C budget is due to the 
dependence of the C yields in the O/H ratio, which 
is higher for Model 4L than for Model 6L. For MS 
the C yields increase with metallicity due to stellar 
winds, while for LIMS the C yields decrease with 
metallicity (Carigi et al. 2005). Since for Models 
6L and 7L the SFR and the O/H values are the 
same, the higher contribution of MS to the C bud- 
get is due to the difference in the IMFs, while for 
Model 6L m up amounts to 80 M Q for Model 7L it 
only amounts to 60 M . Finally the larger contri- 
bution of MS to the C budget in the solar vicinity 
than in Model 1L is due to the evolution of the 
O/H values, and as mentioned before the higher 
O/H values produce higher C yields for MS and 
lower ones for LIMS. 

The difference among the carbon budgets of the 
computed models for NGC 6822 is small and im- 



plies that the result obtained, in the sense that the 
LIMS produce about 63 % of the C abundance and 
that the MS produce about 36 % of the C abun- 
dance, is a robust one. 

7.2. Small infall models 

For the S models we consider the baryonic mass 
aggregation history obtained in section 3.4 un- 
der the assumption of large-scale gas shock heat- 
ing, that prevents a fraction of the gas to be ac- 
creted by the galaxy. Specifically, only 60 % of the 
baryons available after taking into account the ef- 
fect of reionization (see Figures 2 and 3, dotted 
lines) fall into the galaxy. The remaining 40 % of 
the baryons are heated by the collapse of large- 
scale structures (filaments and pancakes) and are 
never accreted by the galaxy. 

In Table 5 we present the main characteristics 
of eight S models. Columns 2 and 3 describe the 
type of outflow and its duration. Columns 4 to 8 
present five properties predicted by the models. 

In Model IS we assume that the gas that falls 
into the galaxy remains in the galaxy during its 
whole evolution. In Figure 11 we show the evolu- 
tion history of six properties of the model up to the 
present time. Model IS fails to fit M gas , since it 
predicts 2.4 times more gas than observed, but re- 
produces all the observational constraints related 
with chemical abundances (see Table 5 and Figure 
11). If the observed M gas value were higher than 
the adopted in this paper, this model would be 
one of the best and would be the only successful 
model without outflow. 

In order to reduce the difference between the 
predicted and observed M gas presented by Model 
IS we have computed models that include well- 
mixed outflows during the early history of the 
galaxy. 

Model 2S is built to reproduce the observed 
Mg as under the assumption of an instantaneous 
well mixed outflow. This event had to occur at 3.1 
Gyr, because at that time the M gas of the galaxy 
amounts to the difference between the observed 
and the predicted M gas by Model IS (2.73 xlO 8 
M ). This model reproduces M gas but the O/H 
value is higher than observed by 3.7 a (see Ta- 
ble 5). The C/H, N/H and Fe/H ratios arc not 
reproduced. 

Model 3S is built to reproduce the O/H ob- 
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Table 4 

Carbon budget for large infall models 1 



Contribution (per cent) 
Model MS LIMS SNIa 



1L 36.3 61.4 2.3 

4L 38.2 59.5 2.3 

6L 37.7 60.1 2.3 

7L 35.9 63.2 0.9 

Solar Vicinity 48.2 49.8 2.0 

Remaining values in the ISM 

1L 36.3 61.4 2.3 

4L 37.2 60.4 2.4 

6L 30.8 67.1 2.1 

7L 34.2 64.8 1.0 

Solar Vicinity 48.2 49.8 2.0 



a Percentage of C in the ISM produced by 
different types of stars over a period of 13 
Gyr. 



served value under the assumption of a well mixed 
outflow that starts at the same time than the star 
formation. To reproduce the observed O/H value 
the outflow phase has to last only 0.1 Gyr. This 
model is very similar to Model IS. During the wind 
phase the outflow rate is equal to 10.0 times the 
SFR. This model reproduces all the observational 
constraints related with chemical abundances, but 
not the M gas observed value (see Table 5). 

Model 4S is built to reproduce M gas value un- 
der the assumption of an early well mixed outflow 
that lasts 1.6 Gyr. Since Model 4S accretes less 
gas than Model 4L, it requires an outflow during 
a shorter amount of time to reproduce the M gas 
observed. During the outflow phase the mass lost 
to the IGM is the maximum possible, that means 
that the gas left is the minimum required to main- 
tain the SFR. 

During the wind phase the outflow rate of 
Model 4S is equal to 8.56 times the SFR. The 
gas mass lost in the outflow amounts to 2.73 xlO 8 
M©, 9.5 times the mass of stars formed during the 
first 1.6 Gyr of the star formation. The predic- 
tions of this model are identical to those of Model 
2S (see Table 5), from ~ 3 Gyr up to now, because 



both models at ~ 3 Gyr had lost practically all of 
their gas and started accumulating gas again at 
the same rate and with the the same SFR during 
their later evolution. In other words, after the end 
of the outflow phase both models have the same 
stellar populations that pollute equally the same 
amount of gas. 

In Model 4S accretion stops after 5 Gyr, when 
the SFR is important. Therefore, its predicted 
O/H value is higher than that obtained by Model 
4L because the chemical elements produced by the 
stars after t= 2.8 Gyr are ejected into an ISM 
which is no longer diluted by primordial infalling 
gas. To lower the O/H value of Model 4S we have 
produced Models 5S - 7S that require even more 
extreme extra assumptions than those needed for 
Models 5L - 7L. 

In Model 5S in order to reproduce the O/H 
value with the least massive cloud we have as- 
sumed that 0.5 Gyr ago a primordial gaseous cloud 
of 1.2 xlO 8 M was captured by NGC 6822. 
With this extra infall, the model can reproduce 
all the observational constraints related to chemi- 
cal abundances with the exception of N/O, but the 
predicted M gas is higher by a factor of 1.6 (see Ta- 
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Fig. 11. — Model IS. Mass formation history and 
predicted chemical evolution. This model assumes 
that 40 % of the gas available after the reonization 
does not fall into the galaxy and this model does 
not include outflows. The present-day observa- 
tional constraints are also presented for compari- 
son. The model fails to reproduce M gas . 

ble 5). The accretion of a huge gaseous cloud of 
primordial material recently is more unlikely than 
at earlier times. If we assume an accretion of a pri- 
mordial cloud at t < 13 Gyr, with the same SFR 
(important at recent times), the model will require 
a more massive cloud producing a larger disagree- 
ment with the observed M gas . We consider Model 
5S more unlikely than Model 5L because the mass 
of the infalling cloud has to be larger and the dis- 
agreement with the observed M gas becomes also 
larger. 

For Model 6S we have assumed a second outflow 
of the selective type during the last 2.5 Gyr, when 
the SFR is most important. In this outflow, 100 % 
of the material ejected by all SNII and SNIa is lost 
to the IGM. This model predicts higher C/O and 
N/O values (see Table 5) because the galaxy keeps 
the C and N produced by LIMS. Model 6S is un- 
likely because between 11.0 and 13.5 Gyr the SFR 
does not indicate the presence of a major starburst 
(see Figure 4) needed to drive the outflow, and the 
potential gravitational well is largest. Moreover, 



Model 6S is more unlikely than Model 6L because 
it requires a recent rich outflow of longer duration. 

In Model 7S we have reduced m up from 80 M 
to 47 M Q and consequently the predicted O/H 
value is reduced by 0.22 dex. This fact produces an 
increase of the C/O and N/O predicted values to 
-0.12 dex and -0.93 dex, respectively, more than 
1.5 a away from the observed values (see Table 
5). The Fe/O value is adjusted because we have 
assumed a lower fraction of binary systems that 
become SNIa (Abin = 0.02). This model does not 
reproduce the C/O and N/O values because with 
this IMF there are more LIMS that are the main 
producers of C and N, at these metallicities. For 
Model 7S the values of M su b, M sta rs, and M rem , 
amount to 3.50 x 10 7 M Q , 1.72 x 10 8 M Q , and 
2.30 x 10 7 M Q , respectively. 
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Fig. 12. — Model 8S. Mass formation history and 
predicted chemical evolution. This model assumes 
the gas mass assembly history of small infall model 
and includes a well mixed outflow during the ini- 
tial 4.0 Gyr of the star formation. The IMF m up 
is set to 60 M Q and the fraction of binary systems 
that becomes SNIa is 2 %. 

No previous S model reproduces all the obser- 
vational constraints. In particular, those models 
that match M gas , predict higher O /H values than 
observed. The extra assumptions included in the 
models produce a rise in the C/O value, higher 
than observed. Consequently, we have computed 
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Table 5 

Inputs and outputs of small infall models 



Model 


Wind 


M gas 




Abundances 






Type 


Duration 


(108 M ) 


12+log(0/H) 


log(C/0) 


log(N/0) 


log(Fe/0) 






(Gyr) 












IS 






4.71 


8.41 


-0.3o 


-1.21 


-1.41 


2S 


Well Mixed 


0.1 


1.99 


8.64 


-0.29 


-1.13 


-1.35 


3S 


Well Mixed 


0.1 


4.51 


8.42 


-0.36 


-1.20 


-1.41 


4S 


Well Mixed 


1.6 


1.98 


8.64 


-0.29 


-1.13 


-1.35 


5S a 


Well Mixed 


1.6 


3.18 


8.48 


-0.30 


-1.14 


-1.36 


6S b 


Well Mixed 


1.6 


1.92 


8.42 


-0.18 


-0.99 


-1.38 


7S c 


Well Mixed 


1.6 


1.98 


8.42 


-0.12 


-0.93 


-1.31 


8S d 


Well Mixed 


4.0 


2.18 


8.50 


-0.22 


-1.04 


-1.43 


Obs. 






1.98 ± 0.20 


8.42 ± 0.06 


-0.31 ± 0.12 


-1.37 ± 0.18 


-1.41 ± 0.10 



a Extra 1.2 x 108 M Q of primordial gas accreted at 13.0 Gyr. 

b All the material expelled by Type II and la SNe during the last 2.5 Gyr is ejected to the IGM. 

c The upper mass limit for the IMF amounts to 47 M Q . Fraction of binary systems that become SNIa amounts 
to 0.02. 

d The upper mass limit for the IMF amounts to 60 M Q . Fraction of binary systems that become SNIa amounts 
to 0.02. 



Model 8S in order to reproduce M gas , O/H, and 
C/O. 

One way to reduce the 0/H and C/H values, 
maintaining the C/0 ratio, is with a well mixed 
outflow of longer duration. In Model 8S, we as- 
sume a well mixed outflow during the first 4 Gyr 
of star formation and the m up of IMF and the 
Abin of Model 7L, the best of the L series models. 
Model 8S reproduces marginally all observational 
constraints with the exception of N/O ( see Table 
5 and Figure 12). 

The C budget for Model 8S is practically the 
same than for Model 7L (see Table 4) because both 
models have the same SFR, IMF, Abin and simi- 
lar O/H. 

During the wind phase the outflow rate of 
Model 8S is equal to 3.13 times the SFR. The 
gas mass lost is 7.2 xlO 7 M Q , 3.5 times the mass 
of stars formed during the first 4.0 Gyr of the star 
formation. Martin (1999) estimates outflow rates 
from 1.7 to 4.8 times the SFR for nearby irreg- 
ular galaxies with the smallest maximum H I ro- 
tation speed of her sample. While these values 
are similar to those required by Model 8S, they 
are smaller than those required by Model 7L. The 
observed outflow rates for outbursting nearby ir- 
regular galaxies seem to support Model 8S over 
Model 7L, but probably at the time of the forma- 
tion of the galaxy the accretion and merging of 
smaller structures also propitiated larger outflow 
rates, it is beyond the scope of this paper to carry 
this discussion further. 

7.3. The N/O ratio 

All our models fail to fit the N/O ratio, this 
could be due to problems with observations and 
with the adopted N yields. In what follows we 
will discuss this problem. 

Carigi et al. (2005) obtained an excellent 
agreement between the present time N/O value 
predicted by the solar vicinity model and the N/O 
values derived from observations of H II regions of 
the solar vicinity. Alternatively, the N/O values 
predicted by all our models are higher than those 
observed in NGC 6822 by Peimbert et al. (2005). 

Moreover the N/O values derived from metal 
poor stars arc in agreement with the model by 
Carigi et al. (2005) but not for metal rich stars, 
these authors suggest that the disagreement is due 



to the N yields adopted in the chemical evolution 
models. 

We consider the M92 yields for C and O for Z Q 
as the best in the literature and we have used them 
in this paper. Unfortunately M92 did not com- 
pute the N yields. In this paper we have adopted 
the N yields for Z Q by MM02. The N yield is 
not self consistent with the C and O yields. Fur- 
thermore since the N abundance is considerably 
smaller than that of C and O, the predictions for 
these two elements are robust, but the N results 
for massive stars should be taken with caution. 

Romano, Tosi, & Matteucci (2005) using N 
yields for LIMS by van der Hoek & Groenewegen 
(1997), different to those adopted by us, cannot re- 
produce the N/O values of NGC 1705, a late type 
dwarf galaxy. On the other hand, Gavilan, Molla, 
& Buell (2005) using their own N yields are able 
to match the present-day N/O value of the solar 
vicinity, but not the N/O evolution. This problem 
should be studied further. 

A possible explanation for the low N/O values 
in H II regions of NGC 6822 could be due to the es- 
cape of an important fraction of the ionizing pho- 
tons. If such were the case, the ionization cor- 
rection factor (ICF) for N becomes considerably 
larger than the usual ICF adopted (Peimbert et 
al. 2005; Relafio et al. 2002). 

8. Discussion 

The chemical evolution of NGC 6822 depends 
on Mtotah the MAH, the baryonic infall as a func- 
tion of time, and the importance of outflows. In 
what follows we will discuss these parameters and 
how they affect the different models that we have 
computed. 

8.1. Total mass 

From the cosmological model adopted we have 
found an M tota i — 2.6 x 10 10 M Q . In what follows 
we will discuss possible upper and lower limits for 
Mtotai- In Figure 13 we present f acc as a function 
of time for the minimum, adopted, and maximum 
Mtotai values. 

The minimum Mtotai f° r the L models is ob- 
tained when there arc no outflows present. For 
the mean MAH and without outflows we obtain 
Mtotai = 1-54 x 10 10 M , for a smaller mass the 
cosmological context predicts a smaller accreted 



21 



0.15 


1 


1 1 I I | 1 1 1 1 | 1 


iii- 
— 

- 


0.1 


- 




_ 

- 






\ 








i\ 












0.05 




V 'M~ al (10'°M o ) = 


2.6 






\ 

\ 


1.54^ '. 
















■ 





1 


i i i i l i i i i l i 


i i i 



5 10 15 

time (Gyr) 

Fig. 13. — f acc , the ratio of the baryonic mass 
accreted by galaxy to the total mass, as a function 
of time for the minimum, the adopted, and the 
high Mtotai values. The triangle, the circle and 
the square represent the ratio of the baryonic mass 
observed to the M to tai for the 1.54 x 10 10 M Q , 
2.6 x 10 10 M Q , and 4.0 x 10 10 M Q , respectively. 

mass than the baryonic mass observed. We have 
computed a chemical evolution model for this to- 
tal mass and find that 12 + log O/H = 8.66 and 
C/O = - 0.29 dex. The O/H value is considerably 
higher than the observed one. To fit the observed 
O /H value we need an O-rich outflow or to reduce 
m up . In both cases the C/O value becomes consid- 
erably larger than observed. For example a model 
with m up of 45 M Q gives 12 + log O/H = 8.42 
in agreement with observations, but it also gives a 
C/O = -0.13 dex, 1.5 a higher than the observed 
value. Therefore the minimum M to tai model can 
be discarded by the chemical evolution models. 
For the S models the minimum Mtotai is smaller 
than 2.6 x 10 10 M Q but larger than 1.54 x 10 10 
M Q . 

The upper limit to M to tai considered by us is 
4 x 10 10 M©. The mass of this model is 1.54 times 
higher than the mass derived in section 3.2. For 
this model the accreted gaseous mass amounts to 
25.1 x 10 s Mq and the ejected mass to 20.8 x 10 8 
M Q . The ejected gaseous mass is about a fac- 



tor of three larger than that of Model 7L and im- 
plies that the efficiency and the duration of the 
outflows have to be considerably larger than for 
Model 7L. If we assume a wind phase of 5.1 Gyr 
duration, like that of Model 7L, the outflow phase 
required amounts to 20.1 times the SFR. This 
value is about an order of magnitude larger than 
that derived by Martin (1999) for nearby irregular 
galaxies. The star formation history, that drives 
the outflow, is observationally fixed, and does not 
change with the increase of the adopted M to tai- 
Moreover the larger M to tai provides a larger grav- 
itational field that reduces the possibility of ex- 
tended outflows. From these arguments we con- 
sider the adopted upper limit for M to tai a generous 
one. 

Based on the models by van den Bosch (2003), 
that predict the virial masses of disk galaxies, 
we can also derive the masses for NGC 6822 for 
the cases with and without outflows. We have 
two observational constraints to derive M to tai, the 
maximum rotational velocity, that amounts to 55 
km s _1 , and the K luminosity that amounts to 
log(L tot /L ) = 8.5. The (L tot /L ) K ratio is ob- 
tained from the NGC 6822 and the solar Mk mag- 
nitudes that amount to -17.91 and 3.33 respec- 
tively (see Table 2, Campins et al. 1985, and 
Hayes 1985). From the two observational con- 
straints and the models by van den Bosch with 
Z = 0.007 and no outflows we obtain M to t al r ^ 
1 x 10 10 M Q , while for the models with Z = 0.007 
and outflow we obtain M tot ai = (2.0 ± 0.5) x 10 10 
M ; these values are somewhat smaller but in 
good agreement with those derived by us. Valcn- 
zuela et al. (2005) find M total 1.93 x 10 10 M© and 
3.4 x 10 10 Mq for their two favored models of NGC 
6822 in good agreement with with our determina- 
tion. 

From the computed chemical evolution models 
we consider that the M to tai adopted by us is rea- 
sonable. But if Mtotai were very different from 
our adopted value, new chemical evolution models 
must be computed. 

8.2. Mass assembly history 

In Figure 1 we have presented three MAHs, the 
slow, the mean, and the fast. In the models pre- 
sented in this paper we have adopted the mean 
MAH. We will discuss why the other MAHs were 
not considered. In Figure 14 we show the f acc for 
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the L models as a function of time for the three 
MAHs under discussion and under the assumption 
that Mtotai = 2.6 x 10 10 M Q . Due to the effects 
of reionization the amount of gas accreted by the 
galaxy reaches 8.5 x 10 8 M Q , 11.7 x 10 8 M , and 
17.0 x 10 8 M Q for the slow, the mean, and the 
fast MAHs respectively. The difference between 
the accreted mass and the observed baryonic mass, 
that amounts to 4.3 x 10 8 M©, has to be ejected 
by outflows. 

The infall of gas produced by the slow MAH at 
early times is not enough to form the older stars 
observed in NGC 6822. Therefore, based on our 
models, we discard the slow MAH. 
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Fig. 14. — face, the ratio of the baryonic mass ac- 
creted by the galaxy to the total mass, as a func- 
tion of time for the slow, mean, and fast MAHs 
for M tota i = 2.6 x 10 10 M©. The filled circle is the 
ratio of the baryonic mass observed to the total 
mass for the three MAHs. 

From Figures 1 and 14 it follows that we require 
to lose a larger amount of mass due to outflows for 
the fast than for the mean MAH models. There- 
fore the fast MAH models would require an early 
outflow of longer duration than the mean MAH 
models, typically of about 8 Gyr. Furthermore the 
gravitational field becomes larger at earlier times 
for the fast MAH than for the mean MAH models 
making it harder for outflows to occur, and the rel- 
ative large initial gas content combined with the 



observed star formation history also prevent the 
presence of large outflow rates. These considera- 
tions made us discard the fast MAH models. 

8.3. Large and small infall models 

As mentioned in section 3.4 several authors 
have suggested that a significant fraction of 
baryons is not accreted by the galaxies due to 
large-scale shock heating of the gas. To consider 
this possibility we have computed a series of S 
models assuming that the evolution up to t = 5.37 
Gyr is the same as for the L models, but after this 
time accretion stops for the S models. The S mod- 
els accrete only 60% of the mass accreted by the 
L models (see Figures 2 and 3). 

The total mass accreted by the S models 
amounts to 7 x 10 8 M©, and since the baryonic 
mass at present is 4.3 x 10 8 M©, these models 
only require to lose 2.7 x 10 8 M©. In Figure 15 
we present f ga i, the ratio of the mass in baryons 
that remains in the galaxy to the total mass as a 
function of time. 

In general the S models predict higher enrich- 
ments in O/H and C/H than the L models (see 
Tables 3 and 5). Model 8S, the best model of the 
S series, requires a well mixed outflow with a du- 
ration of 4.0 Gyr at the beginning of the SFR; 
its O/H, C/H, and M gas values are 1.25a, lAo 
and lcr higher than observed. Based on the chem- 
ical evolution models we favor the large-infall se- 
ries over the small-infall one, particularly we prefer 
Model 7L over Model 8S. 

If the effect of large-scale shock heating on gas 
accretion involves a fraction of baryonic matter 
smaller than that taken into account by the S se- 
ries models, it will be possible to build successful 
galactic chemical evolution models, intermediate 
between Models 7L and 8S. 

9. Conclusions 

We have derived the SFR fort NGC 6822 based 
on observations and photometric evolution mod- 
els. Our photometric models based on the adopted 
SFR, the Salpeter IMF, the KTG IMFs, and three 
Z(t) functions (those predicted by our best chem- 
ical evolution models, 7L and 8S, and that by 
Wydcr) predict absolute magnitudes and colors 
that agree with the observed ones at the one a 
level. This result implies that our derived SFR is 
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Fig. 15. — fg a i, the ratio of the mass in baryons 
that remains in the galaxy to the total mass, as a 
function of time. In the upper panel we show two 
large infall models: Model 1L, without outflows, 
and Model 7L our best model that includes an 
early outflow. In the lower panel we show two 
small infall models: Model IS with no outflows 
and Model 8S that includes an early outflow. Note 
that for Models IS and 1L f gat = f acc . Filled 
circles represent the observational data. 

robust and can be used for different applications 
related to NGC 6822. 

Based on the adopted SFR we present 15 chem- 
ical evolution models for NGC 6822. All mod- 
els evolve according to the mass assembly history 
of non baryonic mass predicted by a cosmological 
context. To distinguish among them we have used 
five additional observational constraints: M gas , 
O/H, C/O, N/O, and Fc/O. 

From the cosmological context adopted and our 
chemical and photometric evolution models for 
NGC 6822 we find that 1.54 < M totai (10 10 M Q ) < 
4.0. The adopted mass for the 15 models presented 
in Tables 3 and 5 is 2.6 x 10 10 M Q . Based on sim- 
ilar arguments, of the three MAHs considered, we 
adopt the mean MAH for all our models. 

For the models of the large infall series, we 
have assumed that during accretion the univer- 
sal baryon fraction is only reduced by reioniza- 



tion. For the models of the small infall series, we 
have assumed that accretion of baryonic matter is 
further reduced by large-scale gas shock heating. 
The best models of each series, Models 7L and 
8S, produce a reasonable fit to the observational 
constraints, but Model 7L is somewhat better. 

Models without well mixed outflows are not 
successful in reproducing the observational con- 
straints because the amount of baryonic gas ac- 
creted by the galaxy is higher than the amount of 
baryonic matter derived from observations. More- 
over based on our adopted SFR the thermal heat- 
ing produced by SNe, under the assumption of no 
outflows, is considerably larger than the binding 
energy of the observed gas content. These two re- 
sults imply that the presence of outflows of well 
mixed material is needed. We argue that these 
outflows most likely occurred early on in the his- 
tory of the galaxy. 

In order to comply with the presence of an early 
outflow, and the observed O/H and C/H values, 
we have computed our best models (Model 7L and 
8S) that have a KTG IMF with m up of 60 M 
while for the solar vicinity m up amounts to 80 Mq . 
That is the number of SN of Type II per generation 
of stars has to be smaller in NGC 6822 than in the 
solar vicinity. A lower value of m up for NGC 6822 
than for the solar vicinity is in agreement with 
recent results on dwarf galaxies by others authors. 

According to Models 7L and 8S and the ob- 
served Fe/O value the fraction of binary systems 
that become SNIa in NGC 6822 is 2 % while in 
the solar vicinity amounts to 5 %. 

The fractions of C produced by MS (m > 8M Q ), 
LIMS (to < 8Mq), and SNIa amount to 60 - 63%, 
38 - 36%, and 1 -2% respectively. On the other 
hand, in the solar vicinity, these values amount to 
48%, 50%, and 2%, respectively. 

The N/O values predicted by all our models 
are higher than observed, the discrepancy could be 
due to the N ionization correction factor adopted 
to derive the N abundance in H II regions and 
to the N stellar yields adopted by our chemical 
evolution models. 

The results obtained in this paper apply to an 
isolated typical irregular galaxy, but do not apply 
to irregular galaxies that have suffered strong tidal 
effects. 

Further progress in modelling NGC 6822 can 
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come if a more precise M tota i is determined, and 
if the issue about the large-scale gas shock heating 
is finally settled. 
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